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Abstract— Pakistan is located just above the tropic of Cancer; 
this offers most optimal locations on the globe for Photovoltaic 
power generation. Khyber Pakhtunkhwa has an average solar 
insolation value more than 5.0KWh/m?/day, which is very 
appropriate for Photovoltaic deployment. But Photovoltaic 
power generation trend shows considerably less progress in this 
region. The aim of this research is to assess the solar potential 
for photovoltaic power generation in Khyber Pakhtunkhwa 
province by identifying feasible sites both technically and 
economically for a utility-scale solar power park installation. 
The feasible sites were identified using Geographic Information 
Systems (GIS) software. This process uses Multi-Criteria 
Analysis method to meet different criteria such as solar 
irradiation, slope and aspect combined with proximity to 
transmission lines and roads plus a number of limiting factors. 
The final results showed that 8% (8,000 km7’) of the research 
area 1s highly suitable for installing utility-scale photovoltaic 
parks. A total of 18 sites with suitability value of 9 and area 
greater than 6 km? have been identified in South and South 
western part of Khyber Pakhtunkhwa. In addition, 70 areas of 
suitability value between 7 and 8 having areas 2-4 km” have been 
identified. Calculations were carried out to find the technical 
potential for power generation. The results showed that 
appropriate amount of feasible areas are available for large scale 
PV installations, with adequate power generation potential. 
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I. INTRODUCTION 


Renewable energies are the fastest growing source of energy 
and considered to be the solution for achieving power system 
sustainability. According to International energy agency 
renewables will provide 30% of the total power demand in 2023, 
as it was 24% in 2017.It 1s forecasted that during the period from 
2018 to 2023, 70% of the world power generation growth will 
be provided by renewables [1]. For solar plant, large area is 
required which causes adverse influence on biodiversity, land 
use, and land cover. Ideally, solar installations should be situated 
on vacant, low productivity agricultural land or on barren land 
and lands covered by bushes to minimize such impact. Non- 


ideal sites are extreme far off location, insecurity and high 
degree of development and covered by forest [2]. The distance 
of solar park from urban population is also very important. The 
Solar plant installed near load centers will have to cover a 
smaller distance and the transmission lines required to transfer 
the energy produced will much less, thus reducing energy loss 
and cost of the energy supplied to consumers [3, 4]. However, 
choosing suitable sites for harvesting solar energy is not a 
straight forward task. It 1s based on multiple sets of physical, 
environmental and socio-economic criteria in order to decide 
optimum geographic locations. Hassan et al. [5] used GIS-AHP 
method for site selection for solar power park in Saudi Arabia. 
Two type of criterion were considered i.e. technical and 
economical. Technical criteria involve are solar irradiation and 
the air temperature. Economic criteria involve slope, aspect, 
Vicinity to power lines, vicinity to cities and vicinity to 
highways. The constraints considered were high slope, land use 
and protected areas comprising of industrial cities, agricultural 
land, national parks, holy places and wildlife sanctuaries. The 
study was divided into four steps. First step was ruling out 
unsuitable sites. In second step an Analytical hierarchy process 
method was used to find out the priority weight and relative 
importance of all criterions. The next step involved overlaying 
different criteria maps as an input criteria according to their 
relative weights obtained in step 2 (AHP) to create an integrated 
analysis. The weighted sum overlay accepts the scaled inputs, 
and adds them together. In final stage a land suitability index 
was developed to show the potential sites suitability distribution 
for solar PV installation. 


Carolina et.al [2] applied Geographical Information System 
(GIS) and multi criteria method to assess the solar energy 
potential for electricity generation in Europe. Factors 
determining the overall suitability (proximity to Electrical grid 
and transportation network, slope, aspect, solar irradiation and 
Population) and different Constrains were identified. Unsuitable 
area identified were forest, wetlands, water bodies, and land use 
whereas undeveloped and short vegetation areas were 
considered as the suitable location. Slope was identified as the 
main suitability factors for land suitability criteria. (Slope 
between 16 and 30 was considered poor while above 30 were 
restricted), solar irradiation (below 900kWh/m/ fall in Poorest 
regions), Proximity to Road network (5000 m cut-off value), the 
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distance to the existing electrical network (grids) and population 
considered as small as possible. All the individual criterion 
factors maps integrated in single layer using weighted linear 
addition technique in GIS to evaluate the total suitability. As a 
validation exercise comparison between existing solar powers 
parks and suitable sites was carried out. The comparison 
indicates that existing solar power parks matched European 
suitability model overlaid in the study. 


Il. METHODODLOGY 


Khyber Pakhtunkhwa, Pakistan is selected as study area on 
the bases of its ideal conditions for solar installation and large 
area. The total area of the Province is 101,741 square km of mix 
terrain, and has annual average irradiance values greater than 
1760KWh/m? according to the National Renewable energy 
laboratory (NREL) solar database. Steps involve are: 


Calculate the total suitable area; 


e Find out how much renewable energy resource is 
present within the calculated areas; 


e Finally calculate the amount of electrical energy that 
could possibly be generated from that available area, 
based on currently available PV technology [6]. 


A. Land suitability criteria: 


Land suitability is based on the following six criteria and 
number of limiting factors. 


i. Solar Irradiation (kWh/m7) 


The NREL is a leading industry for renewable energies 
projects development and installation. The data provided by 
NREL is most credible and is being referenced in literature of a 
lot of scholars. NREL is US based organization that provides 
free data, and have compiled a number of useful datasets for 
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utility-scale or large distributed renewable energy development. 
The Solar data used in this study was acquired from the National 
energy laboratory. When assessing potential sites for PV 
projects, the most important parameter to consider is solar 
radiation levels [6 - 8]. The solar analyst tool in the ArcGIS 
software was used for solar irradiance analysis and mapping [9]. 
The solar irradiation (GHI) for Khyber Pakhtunkhwa is divided 
in 5 classes. Solar irradiation for the study area is illustrated in 
Figure 1. 
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Figure 1. Mean Annual GHI 
Electricity generation from solar energy could be the most 
appropriate source of energy to cope with the rising electricity 
demand because the pattern of global solar irradiation exactly 
follows that of electricity demand. It varies accordingly with 
seasonal variation [4]. Figure 2 shows Seasonal Variation of 
Solar Irradiation. 
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Figure 2. Seasonal Variation of Solar Irradiation 


ll. Slope 


Level plains or minor steep slopes require less construction 
cost as compared to high slopes areas. So, for economic 
feasibility flat terrains are vital for utility-scale PV parks [10]. 
For slope calculation Spatial Analyst toolbar (slope tool) in 
ArcGIS used. For an entire area a slope raster 1s created. This 
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slope raster enables to get an impression of the steepness of the 
terrain, that can be used the output for further analysis. Results 
for slope are illustrated in Figure 3. 
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SLOPE OF KHYBER PAKHTUNKHWA 
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Figure 3. Slope of study Area 


Lands having Slope value less than 5 ° are highly suitable for 
PV installation whereas 5-10° considered as moderately suitable 
[11]. Figure 3 shows the slope of study area. 


lil. Aspect 


Spatial Analyst toolbar in ArcGIS calculates the maximum 
rate of change in value between two neighboring cell. Aspect 
has a direct influence on solar radiation. In northern hemisphere, 
the slopes facing north are mostly shaded, whereas slopes facing 
south receive more solar radiation because these are tilted 
towards the sun and the earth do not directly shade these areas. 
Figure 4 shows result for aspect. 
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Figure 4. Aspect of study Area 


1V. Elevation 


ASTER DEM (30m) has been use for calculating the aspects 
values. Figure 5 shows Elevation of Khyber Pakhtunkhwa. 


V. Distance to Transmission line and Road Network 


Distance of a PV park from transmission lines 1s important 
because of two main reasons; resistance of conductor 1s directly 
proportional to length of the conductor so electrical losses in 
transmission lines increase with the increase in length. Secondly 
constructing new power lines can raise the overall cost of the 
project [12]. The distance between the source and transmission 
line was considered as one key factor [6-8,13,14,15]. The 
distance functions in Spatial Analyst tool was used to find 
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proximity to electrical network. Figure 6, displays the distances 
from transmission line. 


Elevation Map 





Figure 5. Elevation of KPK 
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Figure 6. Proximity to transmission line 
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Figure 7: Proximity to Roads 


Areas < 250 meter are excluded from the study, to counter 
the decrease in module efficiency and damage caused by social 
activities, such as dust due to traffic or any other reason and 
construction activity [16]. Accessibility to potential sites proves 
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to be an important factor. Road accessibility is necessary for 
construction, installation and maintenance throughout the life of 
solar park. Solar power park proximity to a main road 
considered an economic factor [7, 8,12, 14, 15]. Figure 7, shows 
Proximity to Roads. 


V1. Land use and Land cover 


Land cover and land use were considered as constraints in 
this study. Most common restrictions and Constraints applied for 
Utility-Scale solar site are Parks, Water, Wetlands, Forests, 
Wildlife Areas, Urban lands, cultivated and protected land [7- 
11, 14-16]. These constraints are excluded from the suitability 
model. Gaussian Maximum Likelihood algorithm of supervised 
classification was used to extract different features. The result is 
shown in figure 8. 
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Figure 8: Land use and Land Cover 


B. Analytical hierarchy process 


Analytical hierarchy process (AHP) offers the perfect 
compromise solution for contradictory objectives [17]. AHP 
method has been employed to calculate weight for different 
parameters which give importance of a parameter for decision 
making. Table I shows Pair wise comparison matrix. 


TABLE I. PAIR WISE COMPARISON MATRIX 


Elevation | Road Slope | Aspect | DNI Pair wise 
compression 
matrix 
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On applying AHP method [17], weight of each criterion is 
obtained and is shown in Table II. 





TABLE II. PERCENTAGE WEIGHT OF EACH FACTOR 


0.4 1 


0.24 


Category 


0.04 


T/L Network 0.075 





Before overlay each criteria layer is reclassified using 
reclassify tool into classes. When we reclassify a raster, a new 
raster 1s created with the reclassified values. The new value is 
based on aranking scheme usually | to 9. In the third stage, final 
result has been modelled by applying weighted sum overlay 
method using the ArcGIS software. Weighted sum overlay 
technique involves overlaying various factors maps (irradiation, 
slope, distances to road and electrical network) according to 
their relative weights obtained from AHP. A suitability index 
has been prepared to categorize the province on the basis of 
potential solar harvesting site based on the weight of 
contributing factors. Restriction factor like "GHI" less than 4 
KWh/m?/day and lands having slope higher greater than 10° 
were excluded. Slope less than 5° are most suitable [7, 13]. The 
suitability index map ranging from 1-9 has been generated. The 
suitability model is based on six model criteria and restriction 
factors. Table III shows topographic and metrological factor 
rating and Table IV economic factors rating. A priority values 
ranging from 1-9 has been assign to each factors and the 
corresponding values by considering various economic, 
environmental as well as technical aspects. Sites considered 
extremely suitable were assigned a rating of 9, whereas sites 
with the least suitability were set a rating of | or restricted. 


TOPOGRAPHIC AND METROLOGICAL FACTOR RATING 
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TABLE IV. ECONOMIC FACTORS RATING Suitable Sites at N 
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Applying the above rating, restrictions and Weighted -——-—_— 
overlay technique. The suitability model is obtained as shown in 
figure 9. 

Figure 10: Final suitable sites 
Suitability Map N 
= \ TABLE V. SUITABILITY INDEX 
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Figure 9: Suitability Map 


@ Not suitable 


C. Site selection Results 


The Constrains, for instance agricultural land, urban areas, @ Least 
protected areas (National park, holy places etc.) and roads were Mareinall 
excluded from the study. The last step in identifying the best Binaly 
suitable sites was to make sure that the potential sites have m Moderately 
sufficient area for installing utility scale PV Park. Figure 10 


indicates the potential sites for the study area. For utility-scale @ Suitable 
solar PV power park areas having area greater than 1 km? and 
suitability index from 8-9 are consider to be that the potential 
sites [6]. 





Figure 11: Suitability Index 
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The result of MCDM-GIS integration showed that 8% (8000 
km7*) of the concern area is highly feasible for installing solar 
Photovoltaic power parks as shown in figure 11. 


D. Electrical potential Assessment Methodology 
Total Solar Insolation in the suitable area [6]. 
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— /day x 365 (1) 
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Irota/year = Area X 


Trota/ year=8x 10?m?x5kwh/m?x365 
Trota/year=14,600x 10? kwh/year 


Total electrical energy generation 


KWH 
ue I Total x 1] (2) 


Year 


KWH 
= 14,600 x 10°kwh/year x 0.15 


Year 





Total generation=2,190 x10? kWh/year 


CONCLUSION 


The method used in this research work successfully 
identified feasible sites (areas) for photovoltaic deployment 
based on different criteria and multiple stages of evaluation. This 
research work has tried to address the problem in two ways, first 
presenting the results in cartographic form, which provide a 
significant amount of information. Secondly through the use of 
tables and graphs that provides a quantitative view of the results, 
creating a relation between values and their respective locations 
in space. The research consists of two main parts, in first step a 
site suitability analysis was conducted and in second, calculation 
for technical potential and power generation capacity was 
carried out. The overlaid results obtained from the analysis 
showed that 8% area is highly suitable and 10% of the study 
areas have moderate suitability levels. The total insolation for 
suitable area is 14,600x10? kwh/year. This insolation can be 
converted into electrical energy using PV panels having 
efficiency of 15% with optimum tilt angle. The electrical energy 
generation after all the necessary calculation equals to be 2,190 
x10? kWh/year. For utility-scale projects areas less than 1km2 
were excluded from the suitable areas to make sure that total 
size of the land is large enough to be considered for a project. 
The results of the research indicated that, substantial amount of 
suitable sites with adequate power generation potential are 
available for a utility-scale Photovoltaic installation. So it can 
be concluded that the study area is very feasible for on grid or 
off grid system PV deployment. Particularly the area far away 
from the national grid off grid solar PV can be very good choice. 


FUTURE WORK 
This research can be extended in future: 


e Asurvey on social and political acceptance will be very 
supportive. 


e Finding out the solar sites and the effect on voltage 
profile on proximate electrical buses, maximum power 
injection limits and transmission line losses. 
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e The proposed methodology is flexible, it can be easily 
modified and additional physical, environmental and 
planning economic can be included in it. 


e Expand the work to evaluate the economic costs, 
benefits, risks and hazards associated with developing 
the selected sites. 


e The proposed methodology can also use for other 
renewable energy system such as biomass, hydropower 
and geothermal. 
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